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SUMMARY 




The weight required in an initial low Earth parking orbit to accomplish an inter- 
planetary round trip or one-way capture mission is exponentially related to the required 
propulsive velocity increment AV. This total mission AV can be reduced by the selec- 
tion of an efficient parking orbit at the destination planet. This report presents (1) 
several new types of parking orbits that yield low AV's, (2) a comparison of various park- 
ing orbits on the basis of mission AV, and (3) the factors that determine when each type 
of parking orbit is applicable. The study assumes Keplerian trajectories and impulsive 
thrusting. The problems of landing and ascent from the planet surface to the parking 
orbit are not considered. (The term "par king orbit" as used here includes the associated 
approach and departure trajectories as well as the parking orbit itself. ) 

Round trips to Mars in 1979-1980 with total trip durations of 300 to 1000 days were 
selected to illustrate the parking orbits investigated. For these trips, reductions in AV 
of up to 30 percent were achieved by using elliptic rather than circular parking orbits. 

Two efficient means of obtaining the required direction of the Earth return trajectory from 
the parking orbit are (1) by arriving at and departing from the parking ellipse at positions 
away from the ellipse periapsis, or (2) by using a parking ellipse that is out of the plane 
of the interplanetary trajectories and then rotating the plane of the parking ellipse about 
the line of apsides. Elliptic parking orbits will give even larger reductions in AV for 
the more massive planets. 

INTRODUCTION 

The spacecraft weight required in an initial low Earth parking orbit to accomplish 
an interplanetary round trip or one-way capture mission is exponentially related to the 
required propulsive velocity increment AV. The possibility of reducing the mission AV 
by the selection of an appropriate destination planet parking orbit is investigated in this 
report. 

As a basis for comparing the mission AV’s of interplanetary trajectories, it is 
common to assume a low circular parking orbit about the destination planet with the park- 




ing orbit in the plane of the interplanetary trajectories (ref. 1). This destination placet 
parking orbit assumption has the advantage of (1) being similar to the parking orbit fre- 
quently assumed for parking, rendezvous, and assembly at Earth, (2) offering frequent 
opportunities for launching onto the Earth return trajectory, (3) giving low altitude obser- 
vations of the planet surface, and (4) giving the lowest possible atmospheric entry veloc- 
ities for a landing when that is planned. 

The assumption of a low circular parking orbit is adequate for comparing interplan- 
etary trajectories, but when considering the mission itself, the question of an appropriate 
parking orbit must be reexamined. There are several destination planet parking orbits 
that can yield lower mission AV’s than the low circular one (ref. 2); for instance, in 
many cases there is a circular orbit at some altitude above the planet that, of all possible 
circular orbits, yields a local minimum in mission AV. A high circular orbit arrived 
at and departed from by means of semiellipses, which are tangent to the high circular 
orbit and that have a low altitude periapsis, also yields a reduced mission AV. In this 
case the lines of apses of the two semiellipses must be oriented to accommodate the 
arrival and departure interplanetary trajectories (refs. 3 and 4). In the special case 
when the arrival and departure trajectories are tangent, an elliptic parking orbit will 
give still lower mission AV’s. 

When a planetary landing mission is planned, the effect of the parking orbit on the 
landing and takeoff vehicle must also be considered. In reference 4 an analysis was made 
of the initial weight in Earth orbit for Mars landing round trips using the parking orbits 
mentioned previously. In these examples, the elliptic parking orbit that gave the mini- 
mum AV also gave the lowest weight in Earth orbit. Generally there is a wide variety 
of possible parking orbits and associated maneuvers that are potentially of interest and 
that have not been previously systematically studied. The purpose of this paper is then 
to (1) present several new types of parking orbits that yield low AV’s, (2) compare the 
various parking orbits on the basis of mission AV, and (3) illustrate the factors deter- 
mining when each type of parking orbit is applicable. 

To illustrate the comparison between the various types of parking orbits, and to 
show the effect of apsis altitude, the parking orbits are applied to round trips to Mars. 
The trips range from 300 to 1000 days in duration and occur in the 1979-1980 time period. 
The problems of landing and ascent from the surface to the parking orbit are not con- 
sidered. 


SYMBOLS 


e eccentricity 

i, j, k unit vectors along x, y, and z axes, respectively 
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N 

R 

r 

V 

T 

V 

V 
v„ 

AV 

AV c,2 
AV c,3 
AV e,2 
AV e, 3 
AV s 

av t 

x,y,z 

a 

a H,2 

a H,3 

/3 

6 


number of orbits completed 

radial distance from center of force to intersection point of ellipse and hyperbola, 
miles 

radial distance measured from center of force, miles unless otherwise stated 
radius of Mars, miles 
destination planet stay time, days 
velocity, miles/sec 

velocity of Mars in its orbit, miles/sec 

hyperbolic excess velocity of spacecraft (occurs at sphere of influence), 
miles/sec 

propulsive velocity increment, miles/sec 

velocity increment required for entry into circular parking orbit about Mars, 
miles/sec 

velocity increment required for departure from circular parking orbit about Mars, 
miles/sec 

velocity increment required for entry into elliptic parking orbit about Mars, 
miles/sec 

velocity increment required for departure from elliptic parking orbit about Mars, 
miles/sec 

velocity increment savings, miles/sec 
total propulsive velocity increment, miles/sec 
rectangular coordinates 

path angle measured with respect to local horizontal, deg 
vehicle heliocentric path angle at Mars arrival, deg 
vehicle heliocentric path angle at Mars departure, deg 
heliocentric planetary travel angle traversed during stay time, deg 
flight path turning due to gravity, deg 

total gravity turning along hyperbolic trajectories from arriving to departing 
sphere of influence, deg 

flight path turning due to gravity from arriving sphere of influence to initial 
periapsis of Mars centered trajectory, deg 
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6g flight path turning due to gravity from final periapsis to departing sphere of' 

influence, deg 

e angle between Mars velocity vector and local horizontal, deg 

tj true anomaly, deg 

9 total required turning angle measured at Mars sphere of influence, deg 

X position angle of elliptic parking orbit, deg 

X 2 position angle of elliptic parking orbit measured between arrival hyperbolic 

asymptote and elliptic line of apsides, deg 

Xg position angle of elliptic parking orbit measured between departure hyperbolic 

asymptote and elliptic line of apsides, deg 

q o 

p gravitational force constant, miles /sec 

a turning associated with the parking orbit, deg 

r period of parking orbit, circle or ellipse 

cp planetocentric hyperbolic half angle, deg 

oj angular rotation of elliptic orbit about its line of apsides, deg 

Subscripts: 

a apoapsis 

at apo-twist parking orbit 

c circular parking orbit 

e elliptic parking orbit 

H heliocentric 

h hyperbolic orbit 

i indicating either 2 or 3 

lc low circular parking 

op off-periapsis parking orbit 

opt optimum 

p periapsis 

pcd posigrade circularize-decircularize parking orbit 

pe parallel elliptic parking orbit 

R radial distance from center of force to intersection of elliptic and hyperbolic 

orbits, miles 
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r * radial distance from center of force to vehicle, miles 
red retrograde circularize-decircularize parking orbit 
si sphere of influence 

T total 

0 Earth 

cf Mars 

°° sphere of influence 

1 depart Earth 

2 arrive Mars 

3 depart Mars 

4 arrive Earth 


METHOD OF ANALYSIS 

The present analysis considers parking orbits for interplanetary round trips, but it 
can also be applied to one-way capture missions. The propulsive requirements associated 
with various parking orbits depend on the interplanetary trajectory, the characteristics 
of the destination planet, and the planet approach trajectories, as well as the parking 
orbit itself. The subject matter considered herein concerns both the approach and de- 
parture hyperbolas as well as the actual parking orbit. These two items are referred 
to generally as parking orbits. 

The present analysis assumes that the interplanetary trajectory is known; that is, the 
heliocentric velocity vectors for arrival at the planet 2 ’ a H 2 ( ve ^ oc ^ niagnitude and 
path angle, respectively) and for departure V H 3 , 3 are known (fig. 1). This kind 

of information can be generated by the method described in reference 1 , for example. 

The destination planet motion £ (velocity magnitude and path angle) is also assumed 
known. This information may be obtained from an ephemeris, such as reference 5. If 
the planetary motion is known, the planet travel angle /3 is determined by the planet stay 
time, which is defined as the elapsed time between the first AV applied near the planet 
approach periapsis and the AV applied at planet departure. 

To study the parking orbit problem, the coordinate system is changed from one 
centered on the Sun to one centered on the planet by taking the vector difference between 
the vehicle and planet motions (fig. 1). The previous interplanetary trajectory data then 
yield the hyperbolic excess velocity vectors at the planet arrival and departure (V 9 
and 3) as well as the planetocentric turning required (0). The values of ^ 3, 
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Figure 1. - Typical stopover round-trip trajectory showing destination planet en- 
counter both in heliocentric and planetocentric coordinates. 


and 9 depend only on the interplanetary trajectory and are independent of the planet ap- 
proach trajectories and parking orbits. They are the boundary conditions for the study 
of parking orbits. 

In analyzing the planetary parking orbits, several simplifying assumptions are made: 
(1) The planet is spherically symmetric; (2) the maneuvers are made impulsively and the 
trajectories are correspondingly Keplerian; and (3) the minimum radius consistent with 
avoiding atmospheric effects is 1. 1 times the planet surface radius. For many cases the 
parking orbit is assumed to be in the unique plane determined by the arrival and departure 
velocity vectors V M ^ and g and the planet. 

The various parking orbits are compared primarily on the basis of the characteristic 
velocity increment AV required to arrive at, maneuver in, and depart from the parking 
orbit while satisfying boundary constraints imposed by the interplanetary trajectories. 

The AV quite directly affects the initial weight that will be required in Earth orbit to 
perform a mission using the corresponding parking orbit, and a low value is desirable. 
The AV of the parking orbit depends on the apses of the parking ellipse (or radius in 
case of a circular parking orbit), which in turn defines the period of the parking orbit, 
and on the manner by which the turning 9 is obtained. 

The required planetocentric turning 9 is generated by the turning due to the planet 
gravity 6 T and by the turning associated with the parking orbit a. The turning due to 
gravity is the sum of the turning along the approach hyperbola 62 (i. e. , from the sphere 
of influence to the periapses) and the turning along the departure hyperbola 63 (fig. 2 ): 

6 t = 6 2 + 5g (1) 
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where 6 ^ is given by 



The subscript i is a general subscript indicating either 2 or 3. (The basic equations 
required for the present analysis may be found in ref. 3, for example. ) The turning that 
must be supplied by the parking orbit is 


O = e - 6 t (3) 

From equation (3), a positive a indicates that the gravity turning is insufficient. Sim- 
ilarly, a negative a indicates that there is an excess of gravity turning. 

The following sections develop the characteristics of a number of parking orbits, all 
of which supply the turning cr. Unless otherwise stated, all the expressions to follow 
were developed assuming a positive a. For those cases where equation (3) yields a 
negative turning, a should be replaced by (2n - |a|). Several of these parking orbit 
maneuvers are believed to be new; the others are presented for the purposes of compari- 
son and completeness. 


Low Circular Parking Orbit 

The simplest kind of parking orbit is a circular one (fig. 3). This parking orbit 
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should be in that unique plane containing the destination planet and the arrival and depar- 
ture velocity vectors to give the lowest required AV. 

For a circular orbit in general, the total velocity change AV C is composed of two 
parts; the AV required to transfer from the hyperbolic approach trajectory to the cir- 
cular orbit at arrival, AV 0 , anc * the AV required to transfer from the circular to the 

c, z 

hyperbolic orbit at departure, AV 

Cj o 

AV C = AV C> 2 + AV C _ 3 (4) 

where 


AV . = Vu • - V 
c,i h,p,i c 


V. • = & + V 2 
h-P * 1 r 



The period of this parking orbit is given by 



( 6 ) 

( 7 ) 

( 8 ) 


The actual time spent in the orbit, which must be the stay time specified by the 
interplanetary trajectory, can be found from the number of complete revolutions N c 
and the additional travel angle cr (fig. 3). The time spent in the parking orbit, assuming 
a positive, can be written as 



When various parking orbits are compared, the circular parking orbit having a 
radius equal to 1. 1 times the planetary radius will be used as a datum. As mentioned 
earlier, this radius is assumed to be the lowest one outside the sensible planetary atmos- 
phere. From this point forward, "low circular orbit" refers to that orbit having a radius 
equal to 1 . 1 planet radii, and the corresponding AV is designated AV £c . 
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The AV of the remaining parking orbits is compared with that required to achieve 
a low circular parking orbit by a AV savings AV , which is defined as the AV to 
achieve a low circular orbit minus that to achieve the parking orbit under consideration: 

4V g = AV (C - AV parMng orMt (10) 

The parking orbit with the highest AV„ will thus yield the lowest mission AV and is 

s 

hence the most desirable on this basis. The AV e for the low circular orbit is, of course, 

s 

zero. 


Optimum Circular Parking Orbit 


Circular parking orbits are possible at any radius above 1. 1 planet radii. Of these 
orbits, the one which yields the lowest AV is called the optimum circular orbit. For a 
given V ■ a single- maneuver (e. g. arrival) optimum circular orbit may be found by 

°°y 1 

differentiating equation (5) with respect to r and setting the result equal to zero. The 
optimum circular radius found in this manner is 


r 


i, opt 



( 11 ) 


When V^ 2 does not V^ 3 , which is generally the case, two single-maneuver 

optimum circular radii are defined. For this situation, the optimum radius considering 

both the arrival and departure maneuver lies between the two single -maneuver optimums 

and is found by a search in this region of radii. 

The period, time, and AV required for the optimum circular parking orbit may be 

computed by equations (4) to (9) with r. _ , substituted for r. The AV_ may then be 

1 , opt s 

calculated by equation ( 10 ). 


Parallel Elliptic Parking Orbit 

The parallel elliptic parking orbit is illustrated in figure 4 and is defined as follows: 
(1) It lies in the plane formed by V x 2 > 3 , and the destination planet; (2) the ellipse 

and hyperbola periapses are at 1. 1 planet radii; (3) the velocity increments are applied 
at the periapses of the ellipse and hyperbola, which are coincident; (4) thrusting is tan- 
gent to the local velocity; and (5) the major axis or line of apsides of the ellipse is par- 


9 



Figure 4. - Parallel elliptic parking orbit. Elliptic periapsis radius, r^p, 
1. 1 planet radii. 


allel at arrival and departure. These 
qualifications all contribute to a low 
AV T for this parking orbit. 

As was the case for the circular 
parking orbit, the total velocity incre- 
ment for the parallel elliptic parking 
orbit is composed of two parts: 


aV pe = iV e, 2 + aV e, 3 < 12 > 


where 


AV . = V, . - V 
e, 1 h, p, 1 e, p 


(13) 


^ + V 2 


1/2 


h,P,i L ‘ '°°,i 

V r h,p 


(14) 


e,p 


2 ^ r e, a 

r e,p^ r e, a + r e,p) 


1/2 


(15) 


If the results given by equation (12) are used, the AV g can be found from equation (10). 
The elliptic apoapsis velocity, defined now for future use, is 


V = y e > P 
e, a e, p 

e, a 


(16) 


The destination planet stay time is the period of the ellipse times the number of com- 
plete elliptic orbits made: 


T = tN = N JT 
e e 


( r e.a + r e,p)* 

2 M 


1/2 


(17) 


The turning obtained with this parking orbit is only that afforded by gravity. This 
maneuver is therefore limited to those special trips where gravity turning for a periapsis 
of 1. 1 planet radii is equal to the required turning 9 or 


a = 0 


(18) 
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In that the angle a is zero, it is 
expected that this parking orbit will 
yield the highest AV„ of all the maneu- 
vers investigated. It will be shown in 
the DISCUSSION that, for the example 
chosen, this orbit does give the greatest 


Parallel Elliptic Parking Orbit 
by Raising Periapsis 

This parking orbit is like the pre- 
ceding orbit in all the five defining fac- 
tors except number two, the periapsis 

radius. In the preceding parking orbits the periapsis radius was defined to be 1. 1 planet 
radii, and this produces a certain gravity turning. In the case considered now, the peri- 
apsis is raised, and this will reduce the turning afforded by the planet gravity as can be 
seen from equation (2) and as illustrated in figure 5. The turning cannot be increased by 
decreasing periapsis radius because 1. 1 is assumed to be the minimum permissible 
value. Thus, this maneuver applies only when the value of cr for a periapsis radius of 
1. 1 is negative; that is, when gravity gives an excess of turning. 

This orbit is defined by selecting r, , which is equal to r , so that S T = 0, or 

n, p P i 

a = 0. With r known, the AV and period of orbit may be calculated as in the pre- 
e, p s 

ceding section. 



Figure 5. - Reduction ot gravity turning by increasing periapsis radius. 


Apo-twist Parking Orbit 

The apo-twist is one of the new maneuvers considered and is shown in figure 6(a). 

As the vehicle approaches the destination planet, it coasts along a hyperbolic trajectory 
from point A to point B. At B an inplane tangential AV is applied putting the vehicle in- 
to an elliptic parking orbit about the planet; upon reaching point C an additional AV, AV U , 
is applied to change the parking orbit from plane II to plane ID. The elements of the 
ellipse remain constant. A final inplane tangential impulse is applied at B and the depar- 
ture is made along the hyperbolic trajectory B to D. This parking orbit is similar to the 
parallel elliptic orbit in all aspects, except that the plane of the ellipse is out of the plane 
of the interplanetary trajectories. Also, the plane of the parking ellipse is rotated or 


11 


* 



(a) Planetary approach, capture, and departure trajectories. 
Figure 6. - Apo-twist parking orbit. 


twisted about the line of apsides by a thrusting AV W at the ellipse apoapsis; hence the 
name, apo-twist. 

The AV for this type of parking orbit is the same as that for the parallel elliptic 

case with the addition of the twist AV, AV . 

cu 

AV at * AV pe + AV o> < 19 > 

where 

AV a, = 2V e ,a si "(j) <*» 


The elliptic apoapsis velocity is given by equation (16), while the period and time 
spent in the parking orbit are given by equation (17). The twist angle w applied at the 
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apoapsis is found in the following way. 

A coordinate system (fig. 6(b)) is 
chosen so that the line of apses coincides 
with the y-axis. Because the velocity 
vectors V g a are perpendicular to the 
line of apsides, the required twist angle 
o> lies in the x, z-plane. For convenience, 
the departing velocity vector 3 is 
then made to lie in the x, y-plane. The 
twist angle w is also the projection on the 
x, z-plane of the angle between the arriving 
and departing velocity vectors V M g 
3. From the figure it is seen that 

V 3 = V °o, 3 (i cos 6 3 + i sin 6 3 ) (21) 




V x 2 = 2 i cos (180 - u>)cos 6g + j sin 6 2 + k cos 


6 2 sin (180 


u>) 


( 22 ) 


But 


V «, 2 • V », 3 = V », 2 V », 3 cos < 180 - "> (23) 

where i, j, k are the unit vectors along the x, y and z axes, respectively. Solving 
for twist angle gives 


OJ = cos 


-1 


^cos 9 + sin 62 sin 6, 
^ cos cos 63 


( 24 ) 


With the AV at of equation ( 19 ) the AV g can be found as before. 

The apo-twist method can give turning angles only in the direction of increased turn- 
ing, and hence applies only for a positive. For required turning angles that are less 
than those supplied by gravity (a negative), this parking orbit is not applicable. For 
a = 0, this case reduces to the parallel elliptic parking orbit. 
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(a) Planetary approach, capture, and departure trajectories. 
Positive turning associated with parking orbit. 



» 



(b) Planetary approach, capture, and departure trajectories. Negative turning 
associated with parking orbit. 



(c) Effect of off-periapsis thrusting on turning. Negative (d) Detailed geometry for analysis of off-periapsis thrusting. Negative 

turning associated with parking orbit. turning associated with parking orbit. 


Figure 7. - Parallel elliptic parking orbit by off-periapsis maneuvers. Elliptic periapsis radius, r e ^ 1.1 planet radii. 


Off-Periapsis Thrusting Parking Orbit 

The off-periapsis thrusting is the second new maneuver investigated for achieving an 
elliptic parking orbit. It is shown in figure 7(a) for a positive and in figure 7(b) for a 
negative. It has all the characteristics of the parallel elliptic orbit, except that the 
thrusting occurs neither tangent to the local velocity vector nor at the periapsis of either 
the ellipse or hyperbola. While the apo-twist maneuver just discussed is applicable for 
cr positive, this maneuver can be used for o either positive or negative. The following 
development is for a negative. A similar procedure was used for cr positive. 
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(e) Determination of optimum parking orbit orienta- 
tion for two typical Mars missions. 

Figure 7. - Concluded. 


The total AV associated with this park- 
ing orbit is the sum of the AV for arrival and 
for departure: 


AV op= AV e,2 + AV e,3 < 25 > 


To begin, a general description of the 
problem is given. First is considered the 
single maneuver of departing from an ellipse 
oriented at the angle Ag with respect to the 
required departure hyperbolic velocity vector 
V M g (fig. 7(c)). In this instance, if the AVg 
is applied at the ellipse periapsis and tangent 
to the local velocity, the resulting turning ob- 
tained is excessive. However, for periapsis 
thrusting, the required turning can be obtained 
by thrusting at some angle other than zero 
with respect to the local elliptic velocity. 

More generally, the thrusting in addition to 
being noncotangential with the ellipse can occur 


at some elliptic true anomaly r\ 


e, 3' 


For each 


value of 77 e g there is a unique hyperbola 

that yields the required turning, and for a given elliptic parking orbit, g and Ag, one 
value of rj e g will give a minimum AV 0 g. 

This preceding discussion is applicable to the arrival as well as to the departure 
maneuver; that is, for a given A 2 and 2 , 7j g 2 may also be optimized in the manner 
just described to obtain a minimum AV g 2 - However, as discussed later, A 2 and Ag 
are interrelated, and to obtain a minimum, AV Qp = AV g 2 + AV e g, the distribution of angle 
between A 2 and Ag must also be optimized. The following paragraphs describe in de- 
tail first how the individual AV g -'s may be calculated and then how the minimum 

AV op = AV e, 2 + AV e, 3 is found ’ 

The AV required to transfer to or from the elliptic orbit (fig. 7(d)) is, by the law of 
cosines, 


AV 


1 = [ V l R, i + V h, R, i - 2V e , R, i V h, R, i cos ( “ e , R, i " a h, R, i>] 


1/2 


(26) 


The solution of this equation for AV . requires finding the hyperbolic and elliptic veloc- 
ities at the intersection radius R. as well as the angle between them, a . - a, , 

X 6^ Av j X ilj xv ^ X 

where a Q ^ . and ^ . are the path angles measured with respect to the local 
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horizontal of the ellipse and hyperbola, respectively. 

First the terms associated with the ellipse are considered (fig. 7(d)). The ellipse 
will be known in terms of the planet gravitational constant p and its apses (or perhaps 
by its periapsis and period from which the apoapsis may be calculated, eq. (17)). The 
velocity and path angle at point R. (the intersection point) may then be calculated in terms 
of the ellipse true anomaly at point by 


and 


where 


e, R, i 


V 2 + 2p / 1 

e, p 


1 1/2 


R e, i r e, p/ 


a ^=005 -i/WejL 

e ’ R ’ X \V _ .R 

' e, R, i e, iy 


1 + e. 


R„ , = r. 


6,1 e ’ P ' 1 + e e COST ^e,iy 


(27) 


(28) 


(29) 


and e g is the eccentricity given by 


e 


e 


e,p 


e, a 
e,p 


- 1 


+ 1 


(30) 


Next are considered the terms associated with the hyperbola (fig. 7(d)). For the 

hyperbola, the hyperbolic excess velocity V • and the planet gravitational constant 

fj. are known and the ellipse orientation angle is specified. It is desired to find 

V K tj , and a . (eq. (26)) in terms of the hyperbola periapses r. . or true 
it, l n, n, i n, p, i 

anomaly The velocity is given by 


V 


h, R, i 



(31) 


where 
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R e,i 


(32) 


so all the terms are known. The path angle a, . is given by 

n, rt, 1 



where 



(33) 


(34) 


The only unknown term in the relations for R . is the periapsis of the hyperbola 
r h,p,i 

The quantity r, . can be found by the solution of two simultaneous equations in- 
n, p, l 

volving it and the true anomaly of the hyperbola .. It is convenient to write these 
relations in terms of a parametric variable, the eccentricity of the hyperbola e^ which 
also involves r 


h ,P, V 


M + r 


%i 


h,p,i °°,i 


(35) 


From figure 7(d) the first of the two simultaneous equations which results from the angles 
relating the ellipse and hyperbola is 


where 


and 


so that 



(36) 

(37) 

(38) 

(39) 
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The second of the two simultaneous equations is a characteristic equation of the hyperbola: 

1 + e. 


R h, i r h, p, i | 


h,i 


U + e h>i cos 


(40) 


where for the present problem i “ R e i ( 32 ))* Tlie simu l taneous solution of 

equations (35), (39), and (40) gives r h . and 77 ^ ^ this now permits the calculation 

of the path angle a h R . by equation (33) and hence AV g . in equation (26). 

For a given ellipse defined by V . and X-, either r, . or 17 , . defines the 

00 > 1 1 Ilj 1 U, l 

required hyperbola. A search over a range of values of 77 . will reveal the value of 

e, 1 

77 . that yields the lowest AV for the single maneuver, 

e , 1 

Thus far only a single maneuver (arrival or departure) has been considered. How- 
ever, two maneuvers are associated with the parking orbit, the arrival and departure 
(fig. 7(a)). From this figure the A at arrival and departure are interrelated by 


Ag = 77 + 9 - X 2 (41) 

Also from figure 7(a), it can be seen that when 77 9 = 0, the arrival maneuver occurs 

at the ellipse periapsis, and when 77 « = 0 , the departure maneuver occurs at the el- 

e, o 

lipse periapsis. A search is made between these two conditions to find the value of A 2 , 
for instance, that yields the maximum AV savings. 

The AV savings against A 2 for two typical trips is shown in figure 7(e). For each 
point on these curves the values of 77 e 2 anc * r 7 e 3 have been optimized. In each case 
the AV savings maximizes for an orientation of the ellipse within the range previously 
suggested. 

It is not obvious how the stay time at the planet for this parking orbit should be de- 
fined to be most consistent with the previous parking orbits. One possibility is to count 
time from the initial to the final passage of the ellipse line of apsides (fig. 7(a)). In this 
case the period of the orbit is approximately that of the ellipse (eq. (17)) because the 
velocity in the orbit is high near the periapsis, and the time spent in this vicinity is gen- 
erally small compared to the period of the ellipse. 


Perpendicular Apo-twist Parking Orbit 

The perpendicular apo-twist is a special orbit characterized by the line of apsides 
of the parking ellipse being perpendicular to the hyperbola excess velocities V TO ^ anc * 
V^ g. As shown in figure 8 , the vehicle approaches the destination planet by coasting 
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Figure 8. - Perpendicular apo-twist parking orbit. 


along the hyperbolic trajectory from point A to point B. At B an inplane nontangentiai AV 

is applied putting the vehicle into an elliptic orbit about the planet; on reaching point C, an 

additional AV, AV^, is applied to change the parking orbit from plane n to plane HI. 

The elements of the ellipse re main constant. A final inplane nontangentiai impulse is 

applied at B and the planetary departure is made along the hyperbolic trajectory B to D. 

The performance of this maneuver can be calculated by a combination of the apo-twist 

and off-periapsis thrusting case. For this special case, X 0 = x„ = 90° and n «, = 

-o t 6 e,i 

%,3 = ° * 

Also for this case, the turning due to the planet gravity occurs in planes perpendic- 
ular to the plane of 9, and hence, a = Q. To achieve the desired turning, the plane of 
the ellipse is twisted about the line of apsides of the ellipse as in the case of the apo- 
twist maneuver. In this special case the twist angle must be the angle 9 . The twist 
AV is thus 


AV co = 2 V e,a sin (f) < 42 > 

The total AV is the sum of the arrival and departure AV’s as calculated from equ- 
ation (26) and the twist AV. The time in this orbit is given by equation (17). 
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Perpendicular Apo-twist Parking Orbit 
with Optimized Periapsis 

The perpendicular apo-twist with an optimized 
e ' 3 periapsis is a variation of the previous parking 
i orbit where the ellipse periapsis is varied 

(raised) to find that value which yields a maximum 
AV savings. 

Posigrade Circularize-Decircula rize Maneuver 


Figure 9. - Posigrade circuiarize-decircuiarize parking The posigrade circularize-decircularize ma- 

orbit. Elliptic periapsis radius, r ep , 1. 1 planet radn. neuV er is frequently used with an elliptic parking 

orbit to achieve a required parking orbit turning. 
This maneuver has all the characteristics of the parallel elliptic orbit except that the 
lines of apsides of the ellipses at arrival and departure are not parallel (fig. 9). The 
rotation of the line of apsides is obtained by circularizing the arrival semiellipse at its 
apoapsis and then decircularizing into a departure semiellipse (of the same eccentricity 
and periapsis as the arrival ellipse) after the required rotation cr has taken place. This 
rotation occurs in the direction of vehicle motion at the ellipse apoapsis and hence the 
qualification posigrade. 

The total AV for this parking orbit is the AV to arrive at and depart from ellipses, 
given by equation (12), plus the AV to circularize and decircularize the ellipse at apo- 
apsis. This additional AV is twice the AV to circularize the ellipse at apoapsis. The 
total AV is thus 

A Vi = iV pe +2 < V c,a- V e,a> <«> 

where the ellipse apoapsis velocity is given by equation (16) and the circular velocity at 

/ u \ 1/2 

ellipse apoapsis is V = / — t- — l . The AV is again found from equation (10). 

" a UJ 


The minimum time required for the completion of this maneuver is 


T= J ^ g - ’P 

1 2/i 
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* 

Hie time spent at the destination 
planet can be increased by merely re- 



Figure 10. - Retrograde circularize-decircularize parking orbit. 
Elliptic periapsis radius, r e p, 1.1 planet radii. 



Figure 11. - Parallel elliptic parking orbit by turning at sphere of influence. 
Elliptic periapsis radius, r e ^ 1. 1 planet radii. 


maining in one of the elliptic orbits or 
the circular orbit. The stay time will 
then be given by 



where N c and N g are the number of 
completed circular and elliptic orbits, 
respectively. The rotation cr is given 
by equation (3). 

Ret rog rade Ci rcu la rize-Deci rcu la rize 
Maneuver 

The retrograde circularize decir- 
cularize maneuver is a variation of the 
preceding parking orbit where the line 
of apsides is rotated counter to the 
elliptic apoapsis motion (fig. 10). The 
total AV for this parking orbit is 

4V rcd = aV pe +2<V c,a +V e,a) 

(46) 

The minimum time required for 
completion of this maneuver and the 
stay time are obtained from equations 
(44) and (45) with the replacement of 
a by (2 ? t - jcr|). 
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Figure 12. - Variation of heliocentric velocities at Mars arrival and de- 
parture with trip duration. Stopover round trips in 1979-1980. 

Heliocentric velocity of Mars at arrival, 2; heliocentric velocity 
of Mars at departure, 3 heliocentric velocity of spacecraft at 
Mars arrival, V H 2 ; heliocentric velocity of spacecraft at Mars de- 
parture, V H 3 . ’ 

Turn at Sphere of Influence 

In all of the previous parking orbits the required turning a has been achieved in 
close proximity to the planet. In the case considered now (fig. 11), the required turning 
o is made at the sphere of influence. The AV to accomplish a turn is proportional to 
the velocity at which the turn is made. At the sphere of influence the velocity with res- 
pect to the planet, along either the approach or departure hyperbola, is the lowest. Also, 
the lesser of the velocities at the sphere of influence should be chosen. The total AV 
for this parking is then that for the parallel elliptic orbit (eq. (12)) plus the AV for turn- 
ing at the sphere of influence: 


AV si = AV pe + 2Y °o sin ^ ( 47 ) 

2 

The AV savings is given by equation (10) for both ± 0 . The period of the orbit is given 
by equation (17). 
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Figure 13. - Variation of heliocentric path angles at Mars 
arrival and departure with trip dur*ion. Stopover 
round trips in 1979-1980. Vehicle heliocentric path 
angle at Mars arrival, ? vehicle heliocentric path 
angle at Mars departure, « Mars heliocentric path 
angle at Mars arrival, ej % Wars heliocentric path 
angle at Mars departure, ej y 


The parking orbits analyzed in the 
preceding section were incorporated in 
typical Mars stopover round trips (1) to 
illustrate the factors which determine 
when each type of parking orbit is applica- 
ble, and (2) to show a typical comparison 
between the AV savings for the several 
parking orbits. The Mars round-trip 
times range from 300 to 1000 days and 
occur in the 1979-1980 time period. The 
stay times at Mars corresponding to the 
various trip times are given in the abscissa 
of figure 12. The interplanetary trajec- 
tories were selected to give a m in imum 
mission AV assuming a low circular 
parking orbit at Mars. These calculations 
were made using the method of reference 1. 



Total trip time, days 

Figure 14. - Variation of planetocentric hyperbolic excess ve- 
locitities at Mars arrival and departure with trip duration. 
Stopover round trips in 1979-1980. Hyperbolic excess velx- 
ity of spacecraft at Mars arrival, ^ hyperbolic excess 
velocity of spacecraft at Mars departure, V^ y 


Boundary Conditions 

The characteristics of the interplane- 
tary trajectories at Mars in heliocentric 
coordinates are given in figures 12 and 13. 
Figure 12 presents the space ship veloci- 
ties at arrival, V H,2 and at departure, 
V H, g, as well as velocities of Mars, 

g. Figure 13 presents the 
corresponding trajectory and orbit path 
angles with respect to the local horizontal. 
Positive angles are measured clockwise 
from the local horizontal. 

The quantities that are the boundary 
conditions for the planetary parking orbits 
are derived from the heliocentric trajec- 
tory data. They are the hyperbolic excess 
velocities (located at the planet sphere of 
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Total trip time, days 

Figure 15. - Variation of required planetoceritric turning, turning due to gravity, and side of 
periapsis passage with trip duration. Stopover round trips in 1979-1980; elliptic periapsis 
radius, r e 1.1 Mars radii. 


influence) at arrival, 2 , and at departure, 3 , and the turning angle 6. These 
data are presented in figures 14 and 15. The velocities (fig. 14) are highest for the short 
trips, near 300 days. The lowest values occur for the 1000-day trip, which is nearly the 
double Hohmann trip. The required turning 9 (fig. 15) varies from plus 140°, to minus 
100°, and then to plus 175° as trip time increases. The positive direction for 9 is de- 
fined in figure 1 (p. 6). 


Regions of Application 

As was mentioned earlier, the required turning e can be achieved (1) by turning 
afforded by gravity during the approach to and departure from the planet, 6rp = 8 2 + 6^, 
and (2) by turning associated with the parking orbit, a = 6 - 6^. The sign and magnitude 
of <j determined whether some of the parking orbits can be used. 

The turning due to gravity 6 T is greatest for approach hyperbolas with the lowest 
trajectory periapses. The minimum permissible passage radius was assumed to be 
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1 . 1 Mars radii. This is also the periapsis of several of the parking orbits. The gravity 
turning for a 1. 1 radii periapsis is also shown in figure 15. Positive values indicate a 
periapsis on the dark side of Mars, while negative values indicate a periapsis on the 
light side (see sketches on fig. 15). The magnitudes of 6^, for the dark and light side 
locations of the periapses are symmetrical about the zero turn line. 

To minimize the AV associated with any parking orbit, the magnitude of the parking 
orbit turning a should be as s mall as possible. This occurs for every case when 0 and 
6 t are of the same sign. The light or dark side location of the periapsis is, therefore, 
dictated by the sign of Q. With the sign of 0 as a criteria, the location of the periapsis 
alternates from the dark side to the light side and back to the dark side of the planet, 
with increasing trip time, as noted above the abscissa of figure 15. 

When 6^, of the same sign as 0 is used, it can be seen from the figure that 
a = 9 - 6rp is zero for trip times of 485 and 610 days, negative for the trips between 
485 and 610 days, and positive for all other trip times. The sign of a is denoted by the 
type line representing Q. 

Having established the sign of a as a function of trip time, the parking orbits that 
are applicable within each region will be reviewed. The parallel elliptic parking orbit is 
applicable only where o = 0 (periapses of the arrival and departure hyperbola are coinci- 
dent), that is, for 485- and 610-day trip times. The apo-twist maneuver is applicable 
only when a is positive (departure hyperbolic periapsis ahead of arrival hyperbolic 
periapsis), that is, all trip times exclusive of the 485- to 610-day region. The parallel 
elliptic with raised periapsis maneuver is applicable when a is negative, which occurs 
for trip durations of 485 to 610 days. 

Unlike the previous parking orbits, the remaining ones can be used for ±a and, 
hence, for trips of any duration. These parking orbits are the parallel elliptic with off- 
periapsis maneuvers, the perpendicular apo-twist, the sphere of influence turning, the 
low circular, the optimum circular, and both the posigrade and retrograde circularize 
decircularize. 


Comparison of Propulsive Velocity Increment Savings 

It is convenient to use the AV e of three of the parking orbits as references against 

s 

which the AV savings of the remaining parking orbits can be compared. 

Reference parking orbits. - The three reference parking orbits are the low circular, 
the posigrade circularize decircularize, and the parallel elliptic. The first one, the low 
circular parking orbit, was previously selected in the METHOD OF ANALYSIS as the 
datum for calculation of the AV saving. The posigrade circularize-decircularize park- 
ing orbit is of interest because it is frequently considered in the study of round-trip 
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Figure 16. - Variation of Mars capture propulsive velocity increment with elliptic 
parking orbit periapsis and apoapsis radii. Hyperbolic excess velocity of 5 miles 
per second. 


missions. Finally, the parallel elliptic orbit, while in practice applicable only for spe- 
cial trips, is chosen because it gives an upper bound in AV savings. 

An examination of the equations for the AV savings for these three parking orbits 
shows that terms containing the hyperbolic excess velocities do not appear. Thus, these 
parking orbits are independent of ^ anc * V ^ g. Then by assuming that a = 0, the 
AV's for these reference parking orbits become completely independent of trip time. 

The kinds of parking orbits for use as references were selected in the preceding 
paragraph. The periapsis and apoapsis radii of the parallel elliptic and circularize- 
decircularize reference parking orbits remain to be determined. Consider first the 
parallel elliptic orbit. To answer the question of what ellipse will yield the lowest value 
of AV for planetary capture, a typical example was studied. A hyperbolic excess veloc- 
ity of 5 miles per second was assumed and the AV to acquire various ellipses about 
Mars was then computed. The results are presented in figure 16. The AV is plotted 
against the ellipse periapsis for several values of apoapsis radius. For a given value of 
apoapsis, the lowest permissible periapsis gives the lowest AV; a value of 1. 1 Mars 
radii is, hence, used. For a given periapsis radius, the highest apoapsis yields the 
lowest AV. An apoapsis radius of 55 Mars radii was somewhat arbitrarily selected to 
give a 10-day-period ellipse. This will allow several departure opportunities within the 
40-day stay time which was assumed for some of the round trips. All maneuvers having 
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Figure 17. - Comparison of velocity savings for previously 
studied parking orbits with that for reference parking 
orbits as function of trip duration. Mars stopover round 
trips in 1979-19801 


an apoapsis use the 55 Mars radii value. 
The posigrade and retrograde circularize- 
decircularize parking orbits also use a 
periapsis and apoapsis of 1. 1 and 55 Mars 
radii, respectively. 

From the example of the elliptic park- 
ing orbit, it is inferred that in general the 
lowest AV results from an ellipse with 
the lowest periapsis and highest apoapsis 
that are permissible. 

The AV savings for the reference 
parking orbits is given in figure 17 as a 
function of trip duration. The low circular 
orbit is at AV = 0, the posigrade 

b 

circularize-decircularize orbit is at 

1. 21 miles per second, and the upper 

bound for parking orbits restricted by 

r = 1. 1 and r = 55, a parallel elliptic 
P * 

orbit, is at 1. 7 miles per second. 

Previously studied parking o rbits. - 
The parking orbits are discussed in two 
groups, "the previously studied parking 


orbits" and the "new parking orbits”. The previously studied parking orbits are the 
optimum circular parking orbit, turning at the sphere of influence, and the parallel ellip- 
tic with raised periapsis. The actual (as distinct from the reference) low circular and 
posigrade and retrograde circularize-decircularize orbits also f all in this category. The 
actual and reference orbits differ because of the angle a that must be turned in the actual 
case (note that ct = 0 was assumed for the reference cases). The time spent in the actual 
low circular orbit differs from the time spent in the reference circular orbit by the time 
to traverse the angle ct. This time is negligibly small (hours) compared with the stay 
time of 40 days or more. The actual and reference low circular orbit performance are 
thus essentially coincident at all trip times. For the case of the circularize-decircularize 
orbits, the rotation ct occurs in a high circular orbit and can require an appreciable 
amount of time. For this reason the actual and reference parking orbit can differ. If 
the apse radii are specified to be the same, then the AV_ will be the same for the two 
cases but the number of departure opportunities will be fewer for the actual circularize- 
decircularize orbit. If the number of departure opportunities is required to be the same, 
then the apoapsis for the actual parking orbit must be reduced, with a consequent reduc- 
tion in AV savings. 
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Figure 18. - Comparison of velocity savings for new parking orbits with 
that for reference parking orbits as function of trip duration. Mars 
stopover round trips in 1979-1980. 

Of the two circularize-decircularize orbits, the posigrade one will always yield a 
greater AV savings when the apse radii are specified, as shown by figure 17. However, 
if only the number of departure opportunities (or periods) for the parking orbits are 
specified, then the retrograde maneuvers can sometimes give the lower AV. 

With the old maneuvers the largest AV savings is afforded by the posigrade 
circularize-decircularize parking orbit (fig. 17) except for trip times within +15 days of 
where the parallel elliptic parking orbit is possible. The posigrade circularize- 
decircularize parking orbit offers only 72 percent of the upper bound value of AV savings 
and requires four propulsive maneuvers. The new parking orbits were sought to improve 
this situation. 

New parking orbits. - The AV savings for the new parking orbits is compared with 
that for the three reference orbits in figure 18. The new orbits are (1) the apo-twist, 

(2) the parallel elliptic with off-periapsis maneuvers, (3) the perpendicular apo-twist and, 
(4) the perpendicular apo-twist with raised periapsis. The first two parking orbits reduce 
to the parallel elliptic case at 485- and 610-day trip durations. For any particular trip, 
the proper choice of one of these two parking orbits will yield values of AV savings that 
are at least 92 percent of the maximum possible. This is significantly better than the 
72 percent of the maximum possible AV savings obtainable with the reference circularize- 
decircularize parking orbit. In addition, while four thrusting maneuvers are required for 
the circularize-decircularize parking orbit, two and three thrusting maneuvers are 
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Apoapsis in Mars radii, r e< g/r^ 

Figure 19. - Variation (4 velocity savings with apoapsis radius for apo-twist and reference 
circularize-decircularize parking orMs. Stopover round trips to Mars in 1979-1980. 



Apoapsis in Mars radii, r e jrj 

Figure 20. - Variation of velocity savings with apoapsis raidus for parallel elliptic parking 
orbit with off-periapsis maneuvers and for reference circularize-decircularize parking 
orbit. Mars stopover round trips for 1979-1980. 
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required for the parallel elliptic with off-periapsis maneuvers and the apo-twist parking 
orbit, respectively. As discussed previously, the apo-twist parking orbit is not applicable 
to trip durations between 485 and 610 days. 

The perpendicular apo-twist orbits appear to be of interest primarily for the fast 
trips, less than about 450 days, or if a near polar parking orbit is desired. 

Effect of ellipse apoapsis on AV savings. - The preceding comparisons were made 
for parking orbits with an apoapsis of 55 Mars radii. The effect of this assumption is 
discussed with the aid of figures 19 and 20 where the AV savings for several parking 
orbits is plotted against the parking orbit apoapsis radius in Mars radii. The 10-day 
period ellipse is noted on the abscissa. 

The curves of AV„ for the apo-twist and reference circularize-decircularize parking 
orbits, plotted on figure 19, all approach the same asymptote as the apoapsis is increased. 
The shaded band gives the AV savings for all the possible apo-twist maneuvers. It is 
bounded on the upper side by the parallel elliptic orbit (one of the kinds of orbits previously 
used as a reference) and on the lower side by the case of 180° of twist. A 420-day dura- 
tion mission is shown within the band. The reference posigrade circularize-decircularize 
orbit is also shown. The AV savings for the circularize-decircularize and apo-twist 
orbits approach each other for very large apoapsis radii, r a /r^ — °°. Also, the AV 
savings for circularize-decircularize orbit can exceed that for the apo-twist orbit for 
very small apoapsis radii, r a /r^ < 8. However, over the broad middle range of apoapsis 
the apo-twist orbit is superior to the circularize-decircularize orbit. It is also interest- 
ing to note that for the 10- day period ellipse the parallel elliptic AV savings is within 
5 percent of the asymptotic value for which the period is infinite. 

The shaded band of figure 20 gives the AV savings for the parallel elliptic with off- 
periapsis maneuvers. This figure applies for the range of trip times between 485 and 
610 days where the apo-twist cannot be used. This band is bounded on the upper side by 
the parallel elliptic orbit and on the lower side by the condition that a is a maximum. 

For the range of trip times considered, a is a maximum at Q = 0 (from a = 1 9 - 6,p | ) or 
for 540 days (see fig. 15, p. 24). The width of this band is zero for a circular parking 
orbit and increases as the orbit becomes more and more elliptical. An important thing 
to keep in mind about this maneuver is that the AV savings is not independent of the 
hyperbolic excess velocities. Thus, the shaded area shown in figure 20 is valid only for 
the particular hyperbolic excess velocities used, that is, those for the 540-day round 
trip. From this example it is seen that the parallel elliptic orbit with off-periapsis ma- 
neuvers is preferable to the circularize-decircularize parking orbit in the range of apo- 
apsis radii generally of interest. 

Effect of parking orbit on mission AV . - While the primary purpose of this report 
was to compare various parking orbits, it is also of interest to consider the effect of the 
parking orbit on the mission total propulsive AV. This is shown in figure 21 as a function 
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Figure 21. - Variation of total mission propulsive velocity increment 
with trip duration and Mars parking orbit. Mars stopover round 
trips in 1979-1980. Full atmospheric braking at Earth rrturn. 


4 

of trip duration for mission profiles that yield a minimum AV. Because full atmos- 

1 


3 

pheric braking at Earth return is assumed, the propulsion requirement is £ AV. 

1 


The reference parking orbits are represented by the solid lines on the figure. The upper 
line (A) gives the total mission AV required if a circular parking orbit at 1. 1 Mars 
radii is used. The middle line (B) shows the reduction in AV obtained by using the 
circularize-decircularize parking orbit with r n /?A = 1. 1 and r /r /■/ = 55. Finally, 

the lowest solid line (C) shows the lower bound in mission AV as given by the reference 
(a = 0) parallel elliptic parking orbit. Note that on this figure the maximum AV_ of 
1. 7 miles per second is simply the vertical displacement between lines A and C. The 
importance of obtaining a high AV„ is evident from the following examples. The maxi- 
mum AV of 1. 7 miles per second, if obtained, represents a reduction in mission AV 
of 10 percent for the 300-day trip, 21 percent for the 500-day trip, and 33 percent for 
the 1000-day trip. 

The dotted line on the figure shows how closely this minimum AV can be approached 
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using the parking orbits analyzed in this report. As was mentioned previously, 92 percent 
or more of the AV reduction from the low circular to the reference parallel elliptic 
orbit is obtained with the parking orbits herein described. 


Additional Comments 

A few remarks pertaining to the assumptions made in the method of analysis are 
appropriate. The assumption of a spherically symmetric destination planet is equivalent 
to ignoring the effects of oblateness. In actuality, the destination planet’s oblateness and 
resultant shift in the parking orbit orbital plane must be considered. Similarly, the 
rotation of the line of apsides must be accounted for when elliptic parking orbits are used. 
However, it is of interest to note that the effect of oblateness over a given stay time is 
less for orbits with long periods (e. g. , 10 days) than for orbits with short periods (a few 
hours, see ref. 6). 

The assumption of an impulsive or instantaneous velocity change implies an infinite 
thrust. Obviously, in a practical system the thrust will be finite. A low acceleration 
level can enhance the turning due to the planet gravity. Low accelerations are also re- 
lated to low engine weights. These two potential benefits of low thrust must be balanced 
against the reduction in propulsive efficiency (so-called gravity losses) associated with 
low accelerations. 

A final assumption concerns the relation between the interplanetary trajectories and 
the parking orbit boundary conditions. The approach adopted here was to calculate the 
minimum AV interplanetary trajectory for any particular trip time and stay time, and 
then to find the best parking orbit for the resulting parking orbit boundary conditions. An 
alternative is to specify the desired boundary conditions (e. g. , 9 = 6 T ) and then to find 
the minimum AV interplanetary trajectory satisfying these conditions. This approach 
was briefly investigated and found generally to be unattractive in terms of the total mis- 
sion AV. 

In addition to the assumptions discussed previously, several possible areas for further 
study deserve comment. The first of these concerns the manipulation of the various park- 
ing orbits themselves. It is possible that for some boundary conditions a combination of 
maneuvers, say turning at the sphere of influence and apo-twisting, will yield a higher 
AV„ than either single maneuver alone. Second, the parking orbits were evaluated by 
comparing their respective AV savings with no attempt to investigate the interrelation 
between the parking orbit and the problem of landing and taking off from the planet sur- 
face. For instance, the highly elliptical orbit can require a larger vehicle to land a 
given payload on the surface and return it to orbit than does a low circular parking orbit. 
Thus, there is a tradeoff between the AV saving accruing to the interplanetary trajectory 
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by going to a highly elliptical orbit and the resultant increase in the AV for the lan ding 
and takeoff system. Also, the relative difficulty of rendezvousing in the various parking 
orbits must be considered. Both of these effects were evaluated in the Mars missions 
studied in reference 4 and for those missions the elliptic parking orbit gave the lowest 
initial vehicle weight in Earth orbit. 

Finally, numerical examples were shown for only one planet - Mars. The more 
massive planets, such as Venus and Jupiter, may be expected to result in even greater 
variation in AV savings with parking orbit type than does Mars. For Jupiter, for example, 
the AV_ could approach a value of 21.0 miles per second. 

S 


CONCLUSIONS 

A survey of planetary parking orbits applicable to interplanetary round trips and to 
one-way capture missions was made to (1) determine what factors determine when each 
type of parking orbit is applicable, and (2) compare the various parking orbits on the 
basis of their effectiveness in reducing the mission propulsive requirements. Numerical 
results for Mars round trips in 1979- 1980 were presented as examples. The following 
results were obtained: 

1. The applicability of some of the parking orbits depends on whether a (the planet- 
ocentric turning required minus the gravity turning obtained) is positive, negative, or 
zero. 

2. The parallel elliptic parking orbit characterized (1) by thrusting at the elliptic 
periapsis tangent to the local velocity and (2) by an ellipse with the lowest permissible 
periapsis and the highest permissible apoapsis gave the lowest mission propulsive ve- 
locity increment AV. This type of parking orbit is possible only when the turning o is 
zero. For the Mars missions examples, this parking orbit can be used only for trip 
durations of 485 and 610 days. For an ellipse with apoapsis and periapsis of 1. 1 and 

55 Mars radii, respectively, the AV reduction compared with a low circular orbit was 
1. 7 miles per second. This AV reduction is 21 percent of the mission AV for the 
500-day round trip and 33 percent for a 1000-day trip, when full atmospheric braking at 
Earth return is assumed. 

3. When a is positive (insufficient gravity turning), the apo-twist parking orbit may 
be used. This parking orbit can yield 92 percent or more of the AV reduction possible 
with a parallel elliptic parking orbit, and for the example Mars round trips is applicable 
for all trip durations except those between 485 and 610 days. 

4. When a is negative (excess of gravity turning), the parallel elliptic orbit with 
off -periapsis thrusting is attractive. For most trips this parking orbit also yields 

92 percent or more of the AV reduction possible with a parallel elliptic parking orbit. 
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For the Mars trips investigated, it is most useful for trip durations between 485 and 
610 days. 

5. Several other parking orbits were studied. Some are independent of both the sign 
and magnitude of a. None gave as low a mission AV as the parking orbits mentioned 
previously. 

6. It may be expected that the choice of parking orbit will have a larger effect on the 
mission AV for the more massive planets than for Mars. 

7. The final choice of parking orbit will depend on the mission objectives and a com- 
plete systems analysis. In general, however, it is found that, at the expense of additional 
engine firings and possibly other operational complexities, the mission AV may be 
significantly reduced from the usually considered case of a low circular parking orbit. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 4, 1965. 
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